arXiv:1509.01775vl [quant-ph] 6 Sep 2015 


New approach for solving master equation of open atomic system 


Yi-Chong, Ren*’|^and Hong-Yi, Fan'’^ 

^ University of Science and Technology of China, Hefei, People’s Republic of China 
^Shanghai Jiao Tong University, Shanghai, People’s Republic of China 
(Dated: September 8, 2015) 

We describe a new approach called Ket-Bra Entangled State (KBES) Method which enables one convert 
master equations into Schrodinger-like equation. In sharply contrast to the super-operator method, the KBES 
method is applicable for any master equation of finite-level system in theory, and the calculation can be com¬ 
pleted by computer. With this method, we obtain the exact dynamic evolution of a radioactivity damped 2-level 
atom in time-dependent external field, and a 3-level atom coupled with bath; Moreover, the master equation of 
A^-qubits Heisenberg chain each qubit coupled with a reservoir is also resolved in Sec.III; Besides, the paper 
briefiy discuss the physical implications of the solution. 
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I. INTRODUCTION 

The theory of open quantum systems which describes the 
dynamic evolution of system coupled with environment is a 
fundamental approach to understand dissipation and decoher¬ 
ence in quantum opticsC^. One of the most important in¬ 
struments of open quantum systems is represented by mas¬ 
ter equation^^, which plays a significant role in fundamen¬ 
tal aspects, like decoherencel^l^, disentanglement^^l^ , quan¬ 
tum dynamicsl^, non-equilibrium thermodynamics^^, and ap¬ 
plied ones like atom emissionIS, quantum transport^, Bose- 
Einstein condensates. Brownian motionl^, laser field^S] . 

Recently one topic has mostly drawn attention is the 
derivation, solution, and systemic theory analysis of mas¬ 
ter equation which is a set complicated differential equa¬ 
tions of reduced density operator. Generally, the master 
equations are solved by means of phase space representa¬ 
tion of density operator which can transform master equa¬ 
tion into c-number equation such as Fokker-Planck equation 
and Langevin equation ^21. There is also some other way 
like super-operatorE2][ni method, which is applicable only for 
that the Lindblad can be expressed in super-operator genera¬ 
tors of Su(l,l), Su( 2)I^I211 and some other Lie algebra. We 
shall briefly review their method and compare it with our new 
method in Sec.II 

In this contribution, a new procedure called Ket-Bra En¬ 
tangled State(KBES) method have been developed!^, which 
allow one to transform conventional master equation into 
Schrodinger-like equation, then the time evolution of den¬ 
sity operator can be deduced; We can obtain the solution of 
Schrodinger-like equation(see Sec.II). This method is appli¬ 
cable for any master equation of finite-level system in the¬ 
ory. With this KBES method, we obtain the exact solutions 
of three different physical models in Sec.III, which include a 
qubit coupled with reservoir in time-dependent external field, 
a qutrit with Bose-bath, and a A^-qubit Heisenberg chains cou¬ 
pled with reservoir. 


II. SUPER-OPERATOR METHOD VERSUS KBES 
METHOD 

In this section, after briefly reviewing the general super¬ 
operator method, we introduce Fan’s method that constructs 
bosonic thermal entangled representation to mutual transform 
general operator between real and fictitious model^l^ ; As a 
generalization and development of previous, we put forward 
KBES method which can convert operator master equation 
into Schrodinger-like equation. 

To explain the super-operator methodC^HUl ^ consider an 
usual master equation 

P = (1) 

whose Lindblad operator is 

^P “ 2 ~ cr^cr~p - pcr^cr~) 

y 

-f {2cr^pcr~ - cr~cr^p - pcr~cr^) . (2) 

If ^ is time independent, the form solution to Eq. ^ is 

p(t) = e^‘p(0). (3) 

p {t) has the explicit expression if ^ consist of super-operator 
generators of Lie algebras, and there is a wide caiss of mas¬ 
ter equation whose ^ can be expressed in terms of super¬ 
operator generators of Su(2) or Su(I,I) Lie algebra. 

Previous literatures have defined the super-operators 

L+p = cr^per-, L.p = cr-per^, 

Lzp = \ (cr^cr~p - pcr~cr^). ^ 

where [L+, L_] = 2L^, [L^, L+] = ±L+ obeying Su(2) Lie alge¬ 
bra. The explicit expression of ^ is 

^ (2n y l)yy nyL_ + {n+ \)yL+ - yL^. (5) 

Thus Eq. Q can be represented as 

p {f) = ^-\(2n+\)yt+nyL_+{n+\)ytL+-ytL^^ ^ 

_ ^x+{t)L+^\nxJt)L^^xPt)L_ 
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where x+ (t ), (0 is given in Ref. (221 [23l. Finally, the Kraus 

operator solution of Eq. is obtained by transforming the 
super operator L+, L_, into gereral operator 

p{t) = YjKip(!S)K], (7) 


Ki is the Kraus operator and satisfy X; kJKi = 7, where 77, is 
given by 


Kx 

K2 

Kz 

Ka 


+i)r/4 


Xz 

-i)r/4 


X+ (0 (T^, 


(0 cr , 


x+ (0 X- (0 cr^cr , 

[/+ 0-+0-“ (- l)] . 


( 8 ) 


Ulteriorly, prof. Fanl22ll211 constructs the thermal entangled 
State |/) = |0,0) by introducing an extra fictitious mode, 

for the Bose operator 


a\I) = a^I), a^\I) = a\I), (9) 


thus, for the super operator L+ that consist of creation and 
annihilation operator 

L-p\I) = apa^\I) = aap\I) => L_ = aa, 

L+p\I) = a^pa\I) = a^a^p\I) => L+ = a^a\ 

The thermal entangled state |/) can transform general operator 
into fictitious mode, with the characteristic it be used to slove 
master equation of Bose system. 

Next, we shall introduce our new method, consider a den¬ 
sity operator p in Hilbert space 

p = 

m,n 

where |m) constitutes any complete orthogonal basis in 
By introducing an extra fictitious mode, we construct the Ket- 
Bran entangled states 


\r]} = Y^\m,m} ( 12 ) 

m 

fot the density operator p 

|P> = p\n) = n). (13) 

m,n 


Ip) and general operator can mutual transform between real 
and fictitious mode by KBES \t]), which enable one to convert 
master equation into Schrodinger-like equation. 

The case of Eq. Q has been solved by KBES method in 
our previous work ™, here we consider the Lindblad equa¬ 
tion, which is most general form of Markovian and time- 
homogeneous master equation describing non-unitary evolu¬ 
tion of the density matrix p, that is 


dp 

dt 


N^-\ 

^ hn,m[LnpL^ 

n,m=l 

2 ^m^«p)]’ 


(15) 


where H is Hamiltonian part, is arbitrary linear opera¬ 
tors in system’s Hilbert space, and hn,m is constant. With the 
KBES \t]l} = transformed into 

Schrodinger-like equation 


J^\p} = -^{{H-H) + ^Un)\p), 
= ^lIp), 


(16) 


where 


■^lin = 


N^-1 

V h 

n,m=l 





(17) 


Schrodinger-like Eq. ( p^ can be solved by two different ap¬ 
proaches. The first is the evolution operator method, assume 
is independent of time t, then 

\pit)} = e^^‘\pm}. (18) 

where the explicit matrix represent of is demanded, then 
as mentioned before, the calculation of need to diagonal¬ 
ize a -order matrix the calculation is complicate but 
can be done by computer. For few simple situations, we can 
also decompose into several exponential form of oper¬ 
ators with Lie algebra (see Eq. just like super operator 
method. 

The second approach is the stationary state method. Eor 
some high-dimensional systems the calculation of may 
be difficult even with computer, thus we calculate the eigen¬ 
states and eigenvalues of 


^L^i) = (19) 


With the defining of I 77 ), for any operatorA^^^ = \m) {n\m^ 

Amn \rf} = \m, n) = \n) {m\ |??) = \rf}. (14) 

Besides, Eq. is valid for any A = ndmn^nin where 
is real, namel^TA \r]) = A^ \r]). Eq. ( p^ show that the gen¬ 
eral operator can also be transformed between real mode and 
fictitious mode by KBES just as Ean’s method. 

Obviously the fictitious mode indeed represent the bra 
vector of system, so we called [q) Ket-Bra Entangled State 
(KBES). Density operator p can be translated into pure state 


Then the solution can be represented as follow: 

\p{i)) = Yfie"‘'\<fi)^ ( 20 ) 

i 

where Ci is constant which can be determined by initially p ( 0 ) 
and characters of density operator. Obviously, all eigenvalues 
Ai < 0 and the eigenstate \(f) whose eigenvalue A = 0 corre¬ 
sponding the final state p (oo), For t ^ only the eigen¬ 
vector term whose eigenvalue is zero left, other term’s co¬ 
efficient shall disappear along with the growth of r. 
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this characteristic can obtain the final density operator 
directly. Furthermore, most methods of Schrodinger equa¬ 
tion previously can also be used to solve master equation with 
KBES. 

Super-operator method is concise and the solution is ap¬ 
plicable for any initial state, however restricted by Lie al¬ 
gebra, it show a narrow applicable range, and even slightly 
changes of master equation may lead to unsolvable effect for 
super operator method. Compared with the both method, our 
KBES method have three merits: 1. The procedure is ap¬ 
plicable for any master equation of finite-level systems in 
theoretical; 2. The process of resolution is most concise 
and can be fully completed by computer; 3. The method 
can convert master equation into Schrodinger-like equa¬ 
tion, which means most methods of Schrodinger equation 
can be used to solve master equation. These ascendances 
shall be proved in the process of solving follow models. 

III. SEVERAL PHYSICAL MODELS 

To further concretely explain our KBES method, we shall 
introduce three different physical models in this section, more¬ 
over each model can’t be solved by general super-operator 
method. Especially the third model, to best of our knowledge 
is still unsolved for large N, whereas all these can be solved 
by KBES method. 

A. Damped 2-level Atom Model 


where A = ojq - Q and 

Ho = 0.0-^ + Y^ojkalak 

k 

H' = AcTj. +/(cr+ +cr_) 

H'j = (23) 

k k 

Through the derivation of master equation in Ref. ^ , we can 
obtain the corresponding master equation of this model as fol¬ 
low 

^ = -i [H\p\ -r {2cr^pcr~ - cr~cr^p - pcr~cr^) 

1) {2cr~pcr^ - cr^cr~p - pcr^cr ~), (24) 

where H' is the entangled Hamiltonian of qubit and external 

field, y is the spontaneous emission rate and n = n (Q, T) is 

the photon number of frequency Q in temperature T. 

With the corresponding KBES {tjq) = |0,0) -r |1,1), we can 
deduce the Schrodinger-like equation of Eq. ( [24| ): 

^Ip) = ^|p>, (25) 

at 

where a = y{n-\- \) /2,yS = ynl2 and 

^ =i [h' - -r a (2cr~d-~ - cr^cr~ - (T^(T") 

-h p {2(T^d-^ - cr~cr^ - . (26) 


Consider an atom in cavity full of external field that is ra- 
diatively damped by its interaction with the various modes of 
bath in thermal equilibrium at temperature T just as Fig. 1. The 



Fig. 1: Qubit with Cavity Damping in External Field 

Hamiltonian is given in the rotating-wave and dipole approx¬ 
imations as follow: 

H = ccJocr^ + ^CcJk^lak + f(cr^e^^^ 

k 

+ ( 21 ) 

k k 


The formal solution is \p{t)) = ^'^^|p(0)), and the explicit 
formulation of ^ in Kronecker product space can be given: 


' -2a if -if 2p ^ 

if -a-p-2iA 0 -if 

-if 0 -a - p + 2iA if 

. 2a -if if -2p, 


(27) 


where cr- = cr- ® I, = / (g) . In fact, can be obtained 

through the diagonalization of Eq. however the explicit 
expression of p ft) is too long and complicated to be presented 
in here. As we have mentioned in Sec.II, the final density 
p (oo) is represented by the eigenvector of ^ whose corre¬ 
sponding eigenvalue equal to zero. The eigenvector whose 
eigenvalue equals to zero is 


\^) = 


/ (A+/3(Qr+/3))+4A2/3/(a+/3) \ 

l{P+2K^)+{a+pf 
-if{a-l5){a+p-2iK) 
2(/2+2A2)+(a+/3)2 
if{a-l5){a+l5-2iK) 
2(/2+2A2)+(a+/3)2 
(A +a{a+(3)) +4 A^or / {a+p) 

K 2(/2+2A2)+(a+/3)2 . 


(28) 


where / and Q represent the amplitude and frequency of ex¬ 
ternal field respectively. In the interaction picture 

Hi = {H - Ho) 

= H' + H'„ (22) 


Thus p (oo) is given by 


p(oo) 


1 

2^ -h 1 


w H_II_ 

M{nJy,Ap 

ify(n+l /2+2iAy) 

M(n.fj.Ay) 


-ify(n+\/2-2iAy) 'v 
M{n,fy,Ay) 
(n+V)M{nJy,Ky)-f^ 

M{njy,hy) J 


(29) 
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where fy = 2f/y, Ay = 2A/y and 

M{n,fy,Ay) = {n+l/2f + 2 (/^ + 2AI). (30) 

Plot the figure of pn (t) verus yt and Ay on different con¬ 
ditions with initial state p(0) = \<p){<pV 1^) = (|0) + |1)) as 

follow: 



(a)fy = 0,n = 0 (h)fy = hn = 0 (c)fy = 2,n = 0 


Fig. 2: pii versus yt and A^. 


Obviously, pn (oo) that is severely influenced by amplitude 
of external field / and frequency A; For a fixed A, pn (oo) 
show a positive correlation with /, namely pn (oo) increases 
along with the augmentation of /; However for a fixed / > 0, 
there is a negative correlation between pn (oo) and |A| which 
lead to a local maximum of pn (oo) at A = 0; While for / = 0, 
the change of A make no difference to pn (0- For any / 0, 

we find larger |A| shall lead to sharp and frequent fluctuations 
on the outiline of pn (0- 

While the figure of non-diagonal element |pio (t) \ is: Eq. 



(a)/y = 0, ^ = 0 (b)/y = 1,^ = 0 (c)/y = 2, u = 0 


Fig. 3: Ipiol versus yt and Ay. 


( 29| ) show that |pio (oo) p take the maximum 

_ _-'r_ 

Max 


IPlO (oo) 1^ 


4(^+1/2)^-1 + 8/2 


(31) 


when 8A^ + 1 = 2(n-\- 1/2)^ + 4/^. For zero temperature 

n = 0, Ipio (oo) I has a maximum 1/2 V2 at 16A^ y^ = 8/^, 
that has been presented in Fig.3. Despite these, for f = 0, 
change of A make no difference and larger |A| lead to more 
frequent fluctuations on the outiline. 

All those can be easily interpreted , / represents the ampli¬ 
tude of external field, so larger / means strong external field 
which lead to the transition of qubit from ground to excited 
state, i.e. the increasing of pn (oo) and decreasing of poo (oo). 
When A = 0 means the qubit and external field have a same 
frequency coq = ^, there is resonance between qubit and ex¬ 
ternal field, so Pn (oo) have a local maximum value in here; 
The lack of space forbids a further detail discussion of the 
off-diagonal elements pio (0 in here. 


B. 3-level V-type Atom Model 

Now we extend the research of qubit to the case of V-type 
three-level qutrit. Interesting example come from the three- 
level V-type atomic system where spontaneous emission may 
take place from two excited levels to the ground state but di¬ 
rect transition between excited levels is forbidden. However 
the indirect coupling between excited states can appear due to 
interaction with the ground state (quantum interference). 



Fig. 4: V-type Damped Qutrit 

Consider a three-level V-type atom as in Fig.4, in which 
the atom has two nondegenerate excited states |1), |2) with 
transition frequencies to ground state |0) given by cl>i,ccJ 2 and 
spontaneous emission rate yi , 72 respectively. The indirecting 
coupling between |1) and |2) can represent as cross damping 
constant 


712 = Pi (32) 


where Pi change from 0 to 1 represents the mutual orientation 
of transition dipole moments. 

The dynamics of such system is given by the master 
equationl^G^ 

^ (2o-oipo-io - {o-inp}) + ^ (2(To2p(T2o - { 0 - 22 ,p}) 
at 2 2 

7X2 

+ (2croipcr2o + 2cro2pcrio - [(Tn + cr2i ,p}) , (33) 

where [A, B] = AB + BA is the anticommutation and (Tmn 
= \m){n\ means the transition from \n) to \m). Ref. MM 
solve analogous equation with the conditions y^ = 0 and 
7i = 72 - We shall solve this master equation without any 
assumption and approximation in present work. Construct the 
corresponding KBES 


| 77 r) = |00) + |ll) + |22), (34) 

Eq. ( [3^ can be rewrite as Schrodinger-like equation 

^lp> = ^Hp>, (35) 

at 

where 

yi ~ ~ 72 ~ 

(2croi(Toi - crii - d-ii) + — (2cro20^02 “ cr22 “ 0 ^ 22 ) 

712 

+ ( 2 croi(T02 + 2 cr02^01 “ rri2 - 0^21 - cr2l - 0^12) • ( 36 ) 

Then the formal solution is: 


\p{t)) = e^^\pm- 


(37) 
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The explicit solution of Eq. ( [TT] ) can be given by diagonal¬ 
izing to get the explicit expression of which can be 
finished by computer with Mathematica. However the explicit 
expression is too complicated to be given here, so we give the 
solution for ji =72 - the five independent 

matrix elements is: 


Pio (0 = 
P20 (0 = 
P21 (0 = 
Pii w = 
P22 (0 = 


1 ^ ^ ^Pjyt _ P20)j, 

e + P20 + (P20 - Pio)] , 


2 

1 

2^ 


2 

1 

2^ 


P2I -Pl2 + (Pl2 +P2l)^^^^^ 

-(Pii + P22) sinh O^/yO 

Pii - P22 + (pii + P22) cosh (/5/yO 
-(Pi2 +P2i)sinh(yS/yO 

P22 “ Pii + (Pii + P22) cosh (^/yO 
-(Pi2 +P2i)sinh(yS/yO 


(38) 


Remark p^„ = p^„ (0) in the equation for simplify, other 
elements can be easily obtained with pij (t) = p*. (t) and 
P00(0 = 1 -P22(0 -Pll (0- 


When the XXZ chain couple with reservoir as Eq. the 
evolution is no longer unitary and this dynamic process can be 
described by the master equation: 

d ^ 

= -i [Hs,p] + yX i'^cTiPc^i -crtc^iP-PC^tcrj ) (45) 

i=l 

Through the previous models Eq. ( |24|33|) ca n be solved by c- 
number or other method, this moderEq?!?^ is still unsolved. 
To solve the master equation, we construct the corresponding 
Ket-Bra Entangled State 

N 

te> = E Z ,Sn,Si,S2,--- Jn) (46) 

i=l 5 , = 0,1 

where § t = S t = 0,1 represents the spin of i-th qubit. Similar 
with Eq. ( p^ , one can found 

pI'/s ) = IjP<p,^\(P, h = lp>> /47^, 

o-jVls) = d-^\ns)^cr^\ns) = d-~\r}s)- ^ ’ 


C. A/^-qubit XXZ Heisenberg model 

Here we consider a A^-qubit anisotropic XYZ Heisenberg 
chain, the Hamiltonian is given by 

N 

Hs=Yj , (39) 

i=\ 

where {s = x, y, z) are the Pauli matrices of the i-th qubit. 
Js (s = X, y, z) are the strengths of the spin interaction. Eor the 
interaction, when = Jy the model can be called XXZ 
chain. Eq. ( [^ is rewriten as 

= Z +o-r<i)+ ’ (40) 

i 

in which J = + Jy. The eigenvector of Hg can be exactly 

solved by Jordan-Wigner transformation^^ . Now, we only 
consider the ”one-particle” eigenvector (i.e. only one qubit in 
spin up) 


Thus Eq. ( [45] ) can be converted into Schrodinger-like equa¬ 
tion: 

^\p) = ^s\p) = {Hf + G)\p), (48) 

dt 

where 

G = 2yY,cr-d--, 

Hf = i{Hs - Hs)-yYj[o-^o-T +(T+(t 7). 

i 

The formal solution of Eq. ( [48] ) is 

Ip (0) = e^^‘\p (0)> = (0)). (50) 

Introduce the two-mode vector of \k) 

^ I X —^ i 2 n{nk-mk ) ^ 

\k,k'} = \k}{k'\\rjs} = ^ (51) 

m,n=l 


1 ^ 


and the inverse transformation is 




<io“> = T7sE»p 


—i2nnk 


\k). 


For the vector \k) and |0®^), the eigenequation is 


Hs\k} = Ek\k), 
HsW^) = £olO®^), 


the eigenvalues are 


£o = NJ„ 

E^ = Jcos[^) + {N-A)J,. 


(41) 


(42) 


(43) 


(44) 


Easy to found 


where 


G|0,0) = 0, Hf|0,0> = 0; 

G|0, k) = 0, Hf\Q, k) = /to,r|0, k)-, 

G|fc,0) = 0, HF\k,0) = hl^\k,0)\ 

G\k, k') = gk,k' 10,0), HfW, k') = hk,k'\k, k')\ 


gk,k' = 2y6Fk ', 

/!o,r = i (Ek - Eo) - y, 
hk,k’ = i {Ek' - Ek) - 2y. 


With Eq. (52(, for |0, k) and \k, 0) we obtain 


\po,k (t)} = = e’'‘>P\0,k) 

\Pk,o(t)} = = e'^y\k,0} 


(52) 


(53) 


(54) 
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while for \k, k') 

oo ^ 

\Pk,k' (0) = k') = J]^- {h" + GH"-^) \k, k'). 


n=0 


= e'"‘*'‘\k, k'}+^ (A‘" - l) |0,0). (55) 

hk,k' ^ ’ 

Consider the initial state is p (0) = |^)(^| i.e. |p (0)) = \k, k) 
according to Eq. ([5^ and Eq. i 


b (0) = e-^^‘\k, k) + {\- | 0 , 0 ). (56) 

Moreover, fix the initial state p(0) = \ip){ip\,\ip) = + 

where a is real and + \h^ = 1. Accodring to Eq. 
( [4^ and Eq. ( [47] ), the solution can be expressed as 

Ip (0)) = + CLbcr\ + ab*d-\ + |0,0), 

|/?P i2n{k'-k) _ 9 ~ 

- 2^^ ^ 1^,^) + ^2^io,0) 


N 


k,k' = l 


j N N 

at? -link .. at? s 

+ —p ^ 1 ^, 0 ) + —— V ^ | 0 , k) (57) 


As what we have done in Eq. ( [54] ) and Eq. ( [55] ), the solution 
is given by 


\p(t)) = (l - b| 2 e- 2 r') |o, 0 ) + — ^ e^^+''S^V+| 0 ®^, 0 ®'^) 


n,k=l 





Pil—j=l ---j=2,5 ■■■• j=3,4 



Pii— j=l --- j=2,4 ■ ■ j=3 

;|AA A.£k AV\:;/v 

IpAI— j=l ---j=2,4 ■ j=3 


5 10 15 20 

(b) iV = 4 

pil—J=1 ---^2,6 .... j=3,5 M 




Uil—j=l ---j=2,6 .■■■ j=3,5 .-. j=4 


10 15 20 



10 15 20 


(c) iV = 5 


(d) iV = 6 


Fig. 5: p-Jj (0 and Ip^^ (t) | when J -2, y - 1/220, a -\b\-\l V2. 


The solution is given by 


Jtt = 2n7T ± - 


(61) 


m,n,k,k' = \ 


Substituting the value of into Eq. (58), the density p at cross 
points is 

p(t^) = likittmitt) I + 1^1" (l - |0®">(0®^|, 

lb (r*)> = be-^‘"\ip^} + ae-*‘b‘« |0®A, (62) 


+ 


ab* 


1^+10®^, 0®^>. 

n,k=l 


(58) 


Trace to other qubit except 7 -qubit, the density of 7 -qubit is 
obtained 


Poi (f) 

pil (0 


_ ably e^^^+hoj,t 


V2 


(59) 


k,k'=l 


Other elements can be given by PqqO) = 1 - pj^ (0 and 
p{o “ Pqi analysis the decoherence evolution of the 

Heisenberg chain, we plot p^^ (0 and Ip^j (0 | as functions of 
the dimensionless time yt 

In Eig.5, the evolution of pj^ (t) and Ip^^ (0 | is no longer 
monotone decreasing but fluctuate with different frequency, 
the frequency increases with the growth of Jjy. Eig.5 show 
that there are some cross points for A = 3,4, while A = 5 ,6 is 
not; Moreover, the cross points only occurs when A < 4, and 
disappear for A > 5. 

Eor the case A = 4, the abscissa of cross points can be given 
by 


p!i (0=pil it) = pij (0= p\^ it). (60) 


where \ip-) is the IT-states for 4 qubit 

|^±) = 1 (11000) + /|0100) - 10010) + /|0001)). (63) 


That means p {t^) consist of two parts, one is all spin down 
states the other is |<A(^^)) a superposition state of In¬ 

states 1 ^-) and 10 '^'^). p {t^) reduced to a pure IT-states 
for /? = 1,7 = 0. Eq. ( [6T] ) exhibit that the frequency of cross 
points fcp = Jin is independent of 7 , while Eq. ([^ show 
that 7 only influence the exponential damping of both func¬ 
tions. All these imply that the occur of cross point and IT-state 
should be attribute to the interact of Heisenberg model, 7 only 
effects on dissipation process. 

In order that a certain state occurs periodically in the sys¬ 
tem, the necessary and unsufficient condition is that the ratio 
of any two frequencies is a rational number. Eq. ( [44] ) show that 
the condition is satisfied for 2,3,4 and 6 qubits. Eor 6 -qubits 


the Pjj (0 evolve decaying periodically with time, however no 
cross point exist namely 6 -qubit IT-states can’t be created by 
this model. All the solution Eq. ( [5^ , IT-statesEq. ( [63] ), and 
entanglement evolution still need further investigation. How¬ 
ever, in present paper our primary purpose is to introduce our 
new Ket-Bra Entangled State method, so we shall research 
those problems in the near future. 
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IV. CONCLUSION 

In this paper we present a new method (Sec. II) that map 
a master equation into a Schrodinger-like equation, so most 
procedure of Schrodinger equation can be used to solve the 
master equation. All master equation of finite dimension sys¬ 
tem can be resolved by this KBES method in theory. To solve 
master equation of A-level system, the calculation of -order 
matrix’s exponent is necessary (see Eq. ([Tt])), whereas this 
tedious matrix operation can be finished by computer. Eor 
other special cases, the way of stationary Schordinger equa¬ 


tion method may simplify the calculation effectively. 

Through this method, we solve the model of a damped 
qubit in time-dependent external field and a qutrit coupled 
to reservoir, then we resolve A-qubit Heisenberg chain each 
coupled with reservoir at zero temperature, and preliminarily 
analyze the dissipation dynamics and decoherence dynamics 
find that IV-states plays an important role in this process (see 
Sec.III.c). All these cases show that KBES method is a gener¬ 
alization and systematization method for solving master equa¬ 
tion, which can greatly oversimplify the resolution of master 
equation. 
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